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ABSTRACT 



Production of charmed and beauty mesons and baryons in hadron-nucleon collisions 
at not very high energy is considered in the framework of the Quark-Gluon String Model. 
We take into account the possible corrections to the model which disappear at asymp- 
totically high energies but could be important at comparatively low energies, essentially 



in the case of heavy flavour production. The case of nuclear targets is also considered 
and the A-dependence of charmed and beauty hadron yields is calculated. Predictions 
for the case of beauty hadron production at y/s = 39 GeV are presented. 
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1. INTRODUCTION 



Heavy flavour production processes at high energies are considered usually in the 
framework of QCD, see e.g. [1-3], which allows one to describe quantitatively the cross 
sections of c, b and t quark production, their dependence on the transverse momentum of 
the heavy quark, the rapidity distributions, etc. However, these approaches do not allow 
to calculate the spectra of different mesons and baryons. In the case of high energies 
there is also the problem of the parton structure function behaviour in the very small 
x region whereas at comparatively small energies the results depend significantly on the 
values of the heavy quark mass and the QCD scale (see e.g. [4,5]). 

Another approach to heavy flavour production processes has been considered in Ref. 
[6] in the framework of the Quark-Gluon String Model (QGSM). This model is a ver- 
sion of the Dual Topological Unitarization (DTU) and describes quite reasonably many 
features of high energy production processes, including the inclusive spectra of different 
secondary hadrons, their multiplicities, KNO-distributions, etc., both in hadron-nucleon 
and hadron-nucleus collisions [7-10]. High energy interactions are considered as proceed- 
ing via the exchange of one or several pomerons and all elastic and inelastic processes 
result from cutting pomerons or between pomerons [11]. The possibility of exchanging 
a different number of pomerons introduces absorptive corrections to the inclusive cross 
sections which are in agreement with the experimental data on production of hadrons 
consisting in light quarks. The inclusive spectrum of heavy flavoured hadrons in QGSM 
was considered earlier in Refs. [6,12,13]. 

In the present paper we discuss the possibility of using QGSM for calculating heavy 
flavour production at not very high energies. If the ratio mh/y/s is not negligibly small 
some corrections can appear, which disappear with increasing energy. We consider several 
variants of such corrections and compare the results with the experimental data on charm 
and beauty production by proton and pion beams. These corrections can also change 
the predicted A-dependence of charm and beauty production. In particular we present 
quantitative predictions for 5-meson and A& production by protons on nucleon and 
tungsten targets at HERA-B [14] energy, which we assume to be equal y/s = 39 GeV 
(i.e., p lab « 800 GeV/c). 

We also present the results on beauty production at HERA-B energy obtained with 
a QGSM-based Monte Carlo code in which each inelastic collision is considered as a 
hard gluon-gluon one, simulated by PYTHIA. The A-dependence of beauty production 
in hadron-nucleus collisions obtained in this case is very different from two of the three 
variants of the low energy corrections considered in the analytical calculations. 

2. INCLUSIVE SPECTRA OF HEAVY FLAVOURED HADRONS IN QGSM 

As mentioned above high energy hadron-nucleon and hadron-nucleus interactions are 
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considered in the QGSM as proceeding via the exchange of one or several pomerons. Each 
pomeron corresponds to a cylindrical diagram so in the case of cutting a pomeron two 
showers of secondaries are produced, see Fig. 1. The inclusive spectrum of secondaries 
is determined by the convolution of diquark, valence and sea quark distributions u(x, n) 
in the incident particles and the fragmentation functions G(z) of quarks and diquarks 
into secondary hadrons. The diquark and quark distribution functions depend on the 
number n of cut pomerons in the considered diagram. In the case of a nucleon target 
the inclusive spectrum of a secondary hadron h has the form [7]: 

oo 

l/a inel da/dx = J2 w n0n( x ) » (!) 

n=l 

where the functions 4>n( x ) determine the contribution of diagrams with n cut pomerons 
and w n is the probability of this process. Here we neglect the contributions of diffraction 
dissociation processes which are comparatively small in the case of charmed and beauty 
hadron production. 
For pp collisions 

tftt?) = f h qq {x + ,n)f h q {x_,n) + f h q {x + ,n)f h qq {x_,n) + 2(n - l)f?(x +: n)/*(s_, n) , (2) 

x± = - [\J 4mf, Js + x 2 ± x] , (3) 

where f qq , f q and f s correspond to the contributions of diquarks, valence and sea quarks 
respectively. They are determined by the convolution of the diquark and quark distri- 
butions with the fragmentation functions, e.g., 

f q (x + ,n) = [ u q (x 1 ,n)G q (x+/x 1 )dx 1 . (4) 

In the case of a meson beam the diquark contributions in Eq. (2) should be changed by 
the contribution of valence antiquarks. The diquark and quark distributions as well as 
the fragmentation functions are determined via their Regge asymptotics with accounting 
for conservation laws [7,9]. 

In present calculations we use quark and diquark distributions in the proton of the 
form [7]: 



u U u(x,n) — C nu x 2 ' 5 (l — x) n 1-5 , (5) 

Uvd {x,n) = C^-^l-x)"- 1 - 5 , (6) 

u u (x,n) = C u x-°- 5 (l-x) n+0 - 5 , (7) 

u d (x,n) = C d x-°- 5 (l-x) n+1 - 5 , (8) 

Uu(x,n) = u-^(x,n) = CuX~°' 5 [(l + 5/2) 

x (l-x) n+0 - 5 (l-x/3)-5/2(l-x) n+1 ] , n>l , (9) 

u s (x,n) = C s x-°- 5 (l -x) n+l , n > 1 , (10) 
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where 5 = 0.2 is the relative probability to find a strange quark in the sea. The factors 
Cj are determined from the normalization condition 

/ uAx, n)dx = 1 . (11) 
Jo 

The fragmentation functions of quarks and diquarks into charmed mesons and baryons 
were changed a little in comparison with Ref. [12] to obtain a better agreement with the 
existing experimental data. We use these functions in the form 

Gf = GT = ao(l - z) x ^°\l + a lZ 2 ) , (12) 

GZ = Gf = Gf = Gf = G? = Gf = a (l - s^WO , (13) 

Guu = Guu — — — — a o(l - zf +x a * (0) , (14) 

GZ = a (l-zr^\l + a 2 z 2 ) , (15) 

GY = a (l-z) 2+x - a ^°\l-z + a 2 z 2 /2) , (16) 

^ = G^ = a i(l-^) 6+A ^ (0) , (17) 

= G^ = a i(l - z ) 2+x ~ a ^ , (18) 

Gt = G^=Gt(l-z) . (19) 

In the case of A c production there are two different contributions [8]. The first one 
corresponds to the central production of a A C A C pair and can be described by the pre- 
vious formulas. The second contribution is connected with the direct fragmentation of 
the initial baryon into A c with conservation of the string junction. To account for this 
possibility we input into Eq. (2), instead of f q (x + , n) and f q (x-,n), two additional items 
f qq 2(x + ,n) and f qq2 (x-,n) multiplied by (1 — £_) _a *( ) and (1 — x + )~ a *^ respectively. 
The values of f qq 2(x + ,n) and f qq2 (x-,n) are determined by the corresponding fragmen- 
tation functions 

Gt 2 = a 02 z 2 (l-z) 1+x ~^ , (20) 
Gt d2 = a 02 z 2 (l-z) x -"^ . (21) 

The probability for a process to have n cutted pomerons was calculated using the 
quasieikonal approximation [7,15]: 

oo 0. n— 1 

Wn = Onl °n > °n = — (1 ~ e" Z ^ — ) , (22) 
n=l UZ k=0 Kl 

2( ' 7 7Te A « , a P = 87r 7 e A « , £ = ln(s/l GeV 2 ) , (23) 



R 2 + a'i 

with parameters 
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A = 0.139 , a' = 0.21 GeV" 2 , 7 PP = 1.77 , 7^ = 1.07 , 
Rl p = 3.18 GeV" 2 , R\ v = 2.48 GeV~ 2 , C pp = 1.5 , C wp 



In the case of secondary production at not very high energy or if the mass of the 
secondary hadron is comparatively large some problems appear in QGSM. Their origin 
is connected with the normalization (to unity) of quark and diquark distributions. The 
probability for any quark or diquark to fragment into a secondary hadron h is also 
normalized to unity, 

£/ 1 G?( Z )G?(0) t fe = l . (24) 
h Jo 

However the lower limit of integration in Eq. (4) is higher than zero. The minimal value 
of X\ is Xi m i n ~ 171^/ s, so the quark having X\ < x lmin cannot fragment into a secondary 
hadron. As a result the sum of the energy of all secondaries becomes slightly smaller 
than the initial energy. The problem is more clear in the case of a quark which exists in 
the initial hadron. For example in the case of Dp collisions at not very high energy the 
multiplicity of secondary charmed hadrons will be smaller than unity. QGSM is based on 
the Regge theory which contains many other corrections of the same order m\/s which 
have been omitted. However even numerically a small nonconservation of energy and 
quantum numbers can produce problems. To avoid them we consider two possibilities 
to correct QGSM. The simplest way is to use the normalization condition 



f 



1 

Ui(x, n)dx = 1 (25) 



instead of Eq. (11). Another possibility is to multiply every /j in Eq. (2) by a factor 

K = jf fi(x, n)dx I £ ft(x + , n)dx (26) 



and 



D- = f_J?(x,n)dx I f_J?{x-,n)dx . (27) 



3. CHARMED AND BEAUTY HADRON PRODUCTION ON NUCLEON 
TARGET WITH ACCOUNT FOR LOW ENERGY CORRECTIONS 

It is well known that perturbative QCD taking into account leading order (~ a 2 s ) and 
next-to-leading order (~ a 2 ) contributions describes quit well the total cross section 
of heavy flavour production on nucleon targets. However there is a serious problem to 
reproduce simultaneously the experimental cross sections of D-mesons [16] and A c 's [17] 
as well as the shapes of their Xir-distributions. In particular, the Parton Fusion Model 
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[18] assumes that the produced charmed quarks fragment into hadrons thus transferring 
them some fraction of their momentum and cannot recombine with valence quarks be- 
cause of the large rapidity gap between charmed and valence quarks. This model predicts 
the same a^-distribution for favoured and unfavoured .D-mesons in disagreement with 
data [16,19]. Taking into account the intrinsic charm contribution [20] allows one [21] to 
describe the a^-spectra of D-mesons as well as the shape of the A c spectrum. However 
the A c spectrum at moderate xp becomes more than two orders of magnitude smaller 
than the experimental data [17]. 

QGSM does not have so large problems to describe the data [16,17] on D and A c 
production [12]. The largest difference is not greater than 2-3 times, which is not very 
bad keeping in mind the normalization uncertainties. The results of three sets of cal- 
culations, without any corrections (QGSMa), with the corrections (25) (QGSMb) and 
with the corrections (26), (27) (QGSMc), are compared with the experimental data on 
the cross sections for total charm and for different charmed hadrons in Table 1. The 
model parameters for these sets are presented in Table 2. Let us note that the value of 
Oty(0) is comparatively well known [6,12,22] whereas the value of a T (0) can be changed 
significantly. For simplicity we use the same parameter values in the cases QGSMa and 
QGSMb, where the numerical difference of the results is not large. 

Let us compare now the results of our calculations with experimental data (Table 1). 
One can see that the energy dependence of charm production cross section in QGSMb 
and especially in QGSMc cases is not so strong as in QGSMa. Let us note also that the 
QGSMc energy dependence is significantly weaker than perturbative QCD predictions if 
one use modern (i.e., singular at small x) gluon distribution and the same scale values 
at different energies. In all variants we can see a reasonable agreement with the total 
cross section of charm production in pp collisions at 400 GeV/c [16] and 800 GeV/c [23]. 
In the case of SppS energy [24], the experimental errors are too large for excluding any 
variants. In the cases of QGSMa and QGSMb the calculated cross section for D° and 
D + production seems to be too small whereas the cross section of D° production is too 
large. The reasons for this large difference in the calculated values of cr(D) and u{D) 
are the large values of the parameters cq and a 2 in Eq.(8) which is connected with our 
wish to obtain a comparatively large value of <r(A c ) at ^fs = 62 GeV. Correspondingly 
the cross section of A c production at 400 GeV/c [16] is also too large. In the case of 
QGSMc the parameter can not be taken significantly smaller that the used value to 
avoid too small a cross section at \fs = 630 GeV. So the values of a\ and a 2 can not be 
taken so large as in the other cases. As a result the D° and D + cross sections are larger 
whereas cr(A c ) and cr(D) are closer to the experimental data. Let us remark also the 
large difference between old [23] and new [25,26] data on D-meson production at 800 
GeV/c. 

The shapes of the inclusive spectra of .D-mesons do not differ significantly with those 
of [12] as one can see in Fig. 2. 

In the case of beauty hadron production we do not have enough data to fix the 
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parameters of the fragmentation functions. So we use the same fragmentation functions 
as for charm production with parameters 6j (which are also presented in Table 2) instead 
of <2j in Eqs. (12)-(21). Namely for B° production we use the same fragmentation 
function as for D + , for B~ as for D° , for 5° as for D~ , for B + as for D° and for A& as 
for A c . 

The existing experimental point on the total beauty production cross section at y/s 
= 630 GeV [27] allow us to fix the parameter b . The values of bi and b 2 as well as 
&oi and bo2 can be estimated from the condition that the cross section of A c production 
should be (at not very high energy) slightly smaller than the difference of antibeauty 
and beauty production cross section (including the case of a nuclear target). Naturally 
the low energy corrections are more important in the case of beauty production and the 
difference in the predicted energy dependence between the three variants is larger, as one 
can see in Table 3. Actually we normalize the beauty production cross section at \/s = 
630 GeV using the experimental point [27] and obtain very different predictions for low 
energy np collisions. The variants QGSMa and QGSMb predict the last cross section to 
be very small compared to the data [28, 29] and the variant QGSMc predicts it to be 
slightly too large. In the case of beauty production in pp collisions at 800 GeV/c [30] 
the data are in agreement with QGSMa and QGSMb and in strong disagreement with 
QGSMc. However in Refs. [28-30] nuclear targets were used and the values presented in 
Table 3 for nucleon targets are the results of a linear extrapolation. So the experimental 
values possibly can be increased if the A-dependence of beauty production is weaker 
than A 1 . 

Our predictions for beauty production in pp interactions at = 39 GeV (HERA-B 
energy) are presented in Table 4. They are in reasonable agreement with the results of 
the similar calculations of Ref. [13], where the cross section for beauty production was 
equal to 0.01-0.1 fib depending on the parameters. 

4. CHARM AND BEAUTY PRODUCTION ON NUCLEAR TARGETS 

The models based on the DTU approach allow one to calculate the inclusive spec- 
tra of secondaries produced on nuclear targets [31-34] without new phenomenological 
parameters. In multiple scattering theory a hadron-nucleus interaction can be consid- 
ered as a superposition of interactions with a different number of separated nucleons. 
In the QGSM we account for the possibility of multipomeron interaction in every NN 
blob [10,34]. One cut pomeron connects a valence quark-diquark pair of the beam nu- 
cleon with a quark-diquark pair of a target nucleon. All other pomerons connect sea 
quark-antiquark pairs of the incident nucleon with valence quark-diquark pairs or sea 
quark-antiquark pairs of the target nucleons, see Fig. 3. To account for that in our 
calculation of quark and diquark distributions in the beam nucleon the value of n in 
Eqs. (5)-(10) is equal to the total number of cutted pomerons in the nucleon-nucleus 
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interaction, and in the same distributions for a target nucleon the corresponding value 
of n is equal to the number of pomerons connected with a given nucleon. It provides the 
asymmetry of the inclusive spectra of secondaries produced on nuclear targets in compar- 
ison with a nucleon-nucleon interaction. We assume that the secondaries are produced 
ndependently in every shower and their formation zone is large in comparison with the 
nuclear radius. The contribution of intranuclear cascade to the spectra of secondaries, 
as well as the Fermi motion of target nucleons, is neglected. 

A detailed description of light flavour hadron yields from nuclear targets in the QGSM 
was presented in Refs. [10,34]. Now we use the same probability of interaction with a 
given number of target nucleons and consider only the cross sections and the spectra of 
beauty hadrons. 

In the case of QGSMa the charm and beauty cross section in pA collisions has a 
A°-dependence with very small a [10]. It is a consequence of the very strong energy 
dependence of heavy flavour production cross section at comparatively small energies. 
The initial energy is divided between several cut pomerons so the effective energy of every 
shower decreases and the corresponding contribution to charm and beauty production 
cross section is small. So really only the interaction with one target nucleon gives the 
main contribution to the production cross section in QGSMa. In the case of QGSMb 
we have approximately the same behaviour because the corrections are comparatively 
small. However in the case of corrections (26), (27) the situation changes and in this 
case (QGSMc) the value of a becomes closer to unity (see Table 5, where we present 
the calculated values of a for the total charm and beauty production cross section in the 
interval A = l-r-208). 

The predictions of QGSMc for the inclusive spectra of 5-mesons and A^s produced 
in pp and pW collisions at y^i = 39 GeV are presented in Fig. 4. As expected, the cross 
section of all 5-mesons decrease with xf faster than that of A&'s. The yields from a 
nucleus target are larger at small xf and smaller at large Xp, as usual. Some difference 
in the nuclear effects for S-mesons and A&'s can be explained by the diffusion of the last 
ones from the nuclear fragmentation region to the central region as a consequence of the 
comparatively small energy. 



5. MONTE CARLO RESULTS 

In this Section we are going to present the results on beauty production obtained 
with a QGSM-based Monte Carlo code (a detailed description can be found in Ref. 



|p5|). This Monte Carlo model describes hadron-hadron, hadron-nucleus and nucleus- 
nucleus collisions on the same footing, at the partonic level. The nuclear parton wave 
function is constructed as a convolution of parton distributions of individual nucleons 
with the distribution of nucleons inside the nucleus. The position of each nucleon is 
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taken to be described by the Woods-Saxon density: 



p( r ) = Po/(l + exp[(r - r )/a]) , (2£ 

with 



r 



1.19A 1/3 + 1.61A" 1/3 fm, a = 0.54 fm. (29) 



To take the Fermi motion of nucleons into account we generate a Fermi momentum 
p for each nucleon uniformly distributed in the range < p < pp, where pp is the 
maximum Fermi nucleon momentum: 

Vf = (37r 2 ) 1/; V /3 (r) , (30) 

with h = 0.197 fm GeV/c. Isotropical angular distribution is assumed in both the 
coordinate and momentum spaces. 

As to the parton distribution for individual nucleons this is taken to be the same as 
for NN collisions. In particular, the distribution in the number of partons in a nucleon, 
which is directly connected with the value of the multipomeron vertices in the reggeon 
theory (f3~6fl), is taken poissonian: 

w N = exp(-g(s))g N (s)/N\ , (31) 

corresponding to the eikonal picture. The mean number of partons in each nucleon, 
g(s) = goS A , is a function of the center of mass energy y/s. We use go = 3.0 and 
A = 0.09. 

The parton distribution in impact parameter (relative to the center of the correspond- 
ing nucleon) is taken to be gaussian, in accordance with the Pomeron picture of strong 
interactions: 

F(6 p ) = (4vrA)- 1 exp(-62/4A) , (32) 

with the radius depending on the initial nucleon energy. For a projectile or target nucleon, 
A = R 2 + a' In yfe, where a' = 0.01 fm 2 and R 2 = 0.15 fm 2 . 

A hadron or nucleus collision is assumed to be the interaction between partons from 
the projectile and target. A parton from the projectile can interact with one from the 
target if they lie in impact parameter space within an area determined by the parton- 
parton cross section, which has been assumed energy independent, <j v = 3.5 mb. 

In this way the number of inelastic collisions is determined and the inelastic cross 
section is calculated (an elastic event is one with no partons close enough in impact 
parameter space). 

Now, opposite to Ref. |35] (where each inelastic collision could be a hard or a soft one 
with a probability w(s)), each inelastic collision is taken to be a hard gluon-gluon one. 



For simulating these gg collisions the PYTHIA program (|37|) is used. Only gg — > gg 



gg — > qq and gg — ► QQ collisions, with q and Q light and heavy flavours respectively, 
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are considered. In our calculations we have used the EHLQ set 1 structure functions 
(HH), a K = 2.0 factor and a pf n = 2.3 GeV/c. 

Using this model we have computed the xf~ and pt-distributions of particles with 
beauty in pp and pW collisions at \fs = 39 GeV (Figs. 5 and 6). The number of 
generated events is 600000 for pp and 400000 for pW collisions. Also in the last column 
of Table 5 we present the Monte Carlo results for the A-dependence of charm and beauty 
cross sections. 



6. CONCLUSIONS 

In this paper we have presented results of QGSM for heavy flavour production at 
different energies. Three variants of QGSM have been described, which try to take into 
account corrections which are important at not very high energies. Besides, results of a 
Monte Carlo code f35j for production of beauty particles at HERA-B energy have been 



shown. 

We can see that the three variants of QGSM describe reasonably most of the existing 
data on charmed hadron production. Some contradictions are connected with the signif- 
icant differences of different sets of experiment data. For example, the cross section of A c 
production measured in Ref. [16] is very small in comparison with the data of Ref. [17]. 
If the A c production cross section is really so small the difference in D and D production 
cross sections could be obtained smaller by decreasing the parameters a\ and a 2 in Eqs. 
(13), (16) and (17). 

In the case of beauty production at y/s = 39 GeV the model cannot describe si- 
multaneously the data obtained with proton and pion beams [28-30]. As one can see 
in these experimental papers, perturbative QCD also can not describe them with the 
same parameters (QCD scale and quark masses). In our case we can describe the proton 
beam data in the cases of QGSMa and QGSMb and the pion beam data are in rea- 
sonable agreement with QGSMc. A possible explanation of such situation can be the 
assumed A-dependence of beauty production, which could be weaker than A 1 . How- 
ever this would be in disagreement with the direct measurement of the A-dependence of 
.D-meson production [26] where a = 1.02 ± 0.03 ± 0.02. 

So the variants QGSMa and QGSMb predict the beauty production cross section in 
pp collisions at \fs = 39 GeV equal to 5 4- 7 nb in agreement with data [30] and in this 
case they cannot describe the data of beauty production by pion beams [28, 29]. The 
variant QGSMc is in reasonable agreement with the data [29], gives slightly too large 
cross sections in comparison with the data [28] and predicts the cross section of beauty 
production in pp collisions at y/s = 39 GeV equal to 0.2-0.3 fxb with an A-dependence 
~ A om . The A-dependence obtained in the Monte Carlo calculations is A ' 89 , closer the 
the variant QGSMc than to QGSMa or QGSMb. 

We are grateful to R. Rueckl for stimulating discussions. We also thank the Xunta 
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Table 1 



Comparison of the experimental cross sections (in fib) of different charmed hadrons 
produced in pp and irp interactions with the results of our QGSM calculations. In the 
case of np collisions, the cross sections are given only in the region xf > (i.e., in the 
outgoing 71 hemisphere). 



Reaction, y/s 


Experiment 


QGSMa 


QGSMb 


QGSMc 


7r~p -> D + /D , 26 GeV 


5.7±1.6 


9.7 


10.7 


11.2 


vr-p _^ ^0/73°, 26 GeV 


10.1±2.2 


6.3 


7.1 


10.3 


pp -> cc, 27 GeV 


14-23 


28 


33 


30 


pp -> D+/D-, 27 GeV 


11.9±1.5 


12 


14.2 


23.9 


pp -> D°/V, 27 GeV 


18.3±2.5 


23.4 


27.4 


26.6 


pp -> D°, 27 GeV 


10.5±1.9 


3.7 


4.4 


10.5 


pp -> D°, 27 GeV 


7.9±1.5 


19.7 


23 


16.1 


pp -> D + , 27 GeV 


5.7±1.1 


3.7 


4.4 


10.5 


pp — > D~, 27 GeV 


6.2±1.1 


8.3 


9.8 


13.4 


pp -> A c , 27 GeV 


< 6.1 (K/D) 


14.6 


14.6 


7.3 




< 15 (A c /A c ) 








pp — > cc, 39 GeV 


29-55 


47 


52 


41 


pp -> D+/D-, 39 GeV 


33±7 


25 


27 


34 


pp -> D°/D°, 39 GeV 


261S 


40 


44 


38 


pp -> 39 GeV 


17.7±0.9±3.4 








pp-+ D°, 39 GeV 


38±3±13 


6.5 


10 


15.5 


pp -> D + , 39 GeV 


38±9±14 


6.5 


10 


15.5 


pp -> A c , 62 GeV 


(40±18) -r- (204±91) 


19 


19 


8.9 


pp — > cc, 630 GeV 


680±560 ± 250±210 


660 


660 


415 
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Table 2 



The values of the parameters used for the calculations of charm and beauty production 
in QGSM. 



Parameter 


QGSMa 


QGSMb 


QGSMc 


aty(O) 


-2 


-2 


-2 


a 


0.024 


0.024 


0.02 


ai 


10 


10 





a 2 


50 


50 


16 


Ooi 


0.011 


0.011 


0.007 


a 02 


0.005 


.005 


0.0025 


ax(0) 


-8 


-8 


-8 


bo 


0.011 


0.011 


0.0055 


bi 


5 


5 


6 


b 2 


25 


25 


40 


boi 


0.005 


0.005 


0.0015 


bo2 


0.0004 


.0004 


0.0018 



Table 3 



Comparison of the experimental cross sections of beauty production in pp and Tip inter- 
actions with the results of our present QGSM calculations. 



Reaction, \fs 


Experiment 


QGSMa 


QGSMb 


QGSMc 


Ti-p -> bb, 31 GeV 


75±31±26 nb 


2.4 nb 


2.8 nb 


125 nb 


ix~p — > bb, 31 GeV , xf > 


47±19±14 nb 


0.8 nb 


1.0 nb 


63 nb 


n-p -> bb, 34 GeV 


33±11±6 nb 


4.7 nb 


5.5 nb 


154 nb 


pp -> 66, 39 GeV 


5.7±1.5 ± 1.3 nb 


4.8 nb 


6.8 nb 


270 nb 


-> 66, 630 GeV 


19.3±7±9 fjb 


21 fib 


21 fjb 


8 fib 
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Table 4 



Predictions for beauty production in pp interactions at y/s = 39 GeV. All the cross 
sections are in nb. 



Reaction 


QGSMa 


QGSMb 


QGSMc 


pp — > bb 


4.8 


6.8 


270 


pp -> B+/B- 


3.9 


5.5 


220 


pp -> B°/W 


2.1 


3 


210 


PP -> A fe 


1.1 


1.1 


80 



Table 5 



Predictions for the A-dependence (er ~ A a ) of charm and beauty production on nuclear 
targets in the interval A = l-± 208 at y/s = 39 GeV. 



Parameter 


QGSMa 


QGSMb 


QGSMc 


Monte Carlo 


a(cc) 


0.74 


0.75 


0.85 


0.87 


a(bl) 


0.56 


0.57 


0.99 


0.89 
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Figure captions 



Fig.l. Cylindrical diagram which corresponds to the one-pomeron exchange contri- 
bution to elastic pp scattering (a). Its cut which determine the contribution to inelastic 
pp cross section (b). The diagram which correspond to the cut of three pomerons (c). 

Fig. 2. Inclusive spectra of all D-mesons produced in pp interactions at 400 GeV/c 
[16] and its description by QGSMa (solid curve), QGSMb (dashed curve) and QGSMc 
(dashed-dotted curve). The dotted curve represent the result of Ref. [12]. 

Fig. 3. Example of a diagram of interaction of the beam nucleon with two target 
nucleons. 

Fig. 4. Spectra of all 5-mesons and A^'s at ^/s = 39 GeV for a proton beam on 
hydrogen (solid curves) and lead (dashed curves) targets using the QGSMc. 

Fig. 5. :r ^-distributions (normalized to unity) of beauty particles in: a) pp and b) 
pW collisions at y/s = 39 GeV. 

Fig. 6. ^-distributions (normalized to unity) of beauty particles in: a) pp and b) 
pW collisions at y/s = 39 GeV. 
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